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Summary of earlier meetings & problem sets with old (pre 1984) &; new numbering. 


1967 Berkeley 

1968 Berkeley 

1969 Asilomar 


1970 Tucson 

1971 Asilomar 

1972 Claremont 

72:01-72:05 

1973 Los Angeles 

73:01-73:16 

1974 Los Angeles 

74:01-74:08 

1975 Asilomar 

75:01-75:23 



1976 San Diego 

1-65 i.e., 

76:01-76:65 


1977 Los Angeles 

101-148 i.e., 

77:01-77:48 


1978 Santa Barbara 

151-187 i.e., 

78:01-78:37 


1979 Asilomar 

201-231 i.e., 

79:01-79:31 


1980 Tucson 

251-268 i.e., 

80:01-80:18 
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301-328 i.e., 

81:01-81:28 


1982 San Diego 

351-375 i.e., 

82:01-82:25 


1983 Asilomar 

401-418 i.e., 

83:01-83:18 


1984 Asilomar 

84:01-84:27 
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1986 Tucson 
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1987 Asilomar 

87:01-87:15 

1988 Las Vegas 

88:01-88:22 

1989 Asilomar 

89:01-89:32 
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1991 Asilomar 

91:01-91:25 

1992 Corvallis 

92:01-92:19 

1993 Asilomar 

93:01-93:32 

1994 San Diego 

94:01-94:27 

1995 Asilomar (present set) 95:01-95:19 


[With comments on: 90:14, 94:07, 94:12, 94:13, 94:15, 94:26] 

UPINT(2 ) = Richard K. Guy, Unsolved Problems in Number Theory, Springer, 1981. (Second 
edition 1994). 

COMMENTS ON ANY PROBLEM WELCOME AT ANY TIME 

Centre for Number Theory Research, 

MPCE Building E7A, 

Macquarie University, 

NSW 2109 Australia. 

gerry@mpce.mq.edu.au Australia-2-9850-8952 fax 9850-8114 

15 Aug 96 


1 









COMMENTS ON EARLIER PROBLEMS 

90:14 (D.H. and Emma Lehmer) Let p = n 2 + 108 = 6/ + 1 = 109, 157, 229, 277, 397, ... be 
a prime, and 

’h=Y.xr 

i =0 

(C P is a primitive j>th root of unity) be the Gaussian periods of degree 6. Let 
Sj = rjj - ij j+k k 16 

(i.e., k is prime to 6) Then the 5 :j satisfy the sextic 

G(x) = x 6 — p(x — l) 2 (x — 2) 2 = 0 

so that 

5 

<2(1) = life “ -hi+fc - 1) = 1 

i=o 

and hence pj = rjj — rjj +k — 1 is a unit. Is pj a fundamental unit? [These primes do not have 
translation units.] 

For e = 5 there is only one prime, p = 211, known to have this property. 

[Editor’s note: I think this means that 211 is the only prime of the form 5/ + 1 for which it 
is known that rjj — r/ J+fe — 1 is a unit, where rjj is a Gaussian period of degree 5.] 

94:07 (Bart Goddard) Let 

fn{x) = (n — \)\x n + (n — 2)\x n ~ 1 +-1- 2x 3 + x 2 — x H-j- 

and f n {a n ) = i»f [0 ,,] f„{x) 

(a) Find an upper bound for f n (a n )- 

(b) Find lim^^ a n 

Remark: Possible answers: (a) < y/n, (b) 

Paul Feit gives the following solution to (b): 

We answer Part (b) by proving that lim^ooCtn = 0. To begin, consider the derivatives of the 
functions: 

f' n = n\x n ~ l +-b 2x — 1. 

Each derivative f n can be expressed as g n {x) — 1, where g n (x) is a polynomial with no constant 
term and only non-negative coefficients. In particular, on [0, oo), g n {x) assumes 0 at x = 0 and 
increases. 

By elementary calculus, it follows that a n can be characterized by the properties 


1. g n (x ) < 1 for 0 < x < a n , and 

2. g n (x ) > 1 for a n < x < 1. 




Clearly a n < 1/2 for all n. To show that litn»^oera« = 0, it suffices to prove that, for any r € (0,1], 
there is an N such that g n (r) > 1 for every n > N. In fact, since g n (X) is a sum of positive 
monomials, it suffices to prove the inequality with a single term in place of g n (x). 

Standard theory proves that for any r > 0, 

lim„^ 00 n!r n_1 = oo 

In the present context, this remark proves Part (b). 

This method will produce estimates for f n {a n ) as well. (That is Part (a) of the problem.) 
The derivative of f n stays between 0 and 1 on (0, a n ). Consequently, f n (a n ) must lie between 
/„(0) = n " and n ” — a n . For n large, the latter is a very narrow range. 

94:12 (Graeme Cohen, via Doug E. Iannucci) For a positive integer n, show that a(a(n)) = 4n 
only if n is odd. 

Remark: The only such numbers less than 5 • 10 8 are 15, 1023 and 29127. Graeme Cohen 
has further checked this to 10 9 . The results will appear in 

G. L. Cohen & H. J. J. te Riele, Iterating the sum-of-divisors function, 
which was presented at CANT’95 at Macquarie University, Sydney, last April. These numbers 
were mentioned in 

Carl Pomerance, On multiply perfect numbers with a special property, Pacific J. Math., 57(1975) 
511-517. 

and Carl Pomerance sends the further references: 

G. G. Dandapat, J. L. Hunsucker & C. Pomerance, Some new results on odd perfect numbers, Pacific 
J. Math., 57(1975) 359-364. 

J. L. Hunsucker & C. Pomerance, There are no odd super perfect numbers less than 7 • 10 24 , Indian 
J. Math. 17(1975) 107-120. 

In the first it is shown that if n is odd and super perfect, then neither n nor <r(n) is a prime 
power, and that either n or a(n) has at least 3 distinct prime factors. In the second the title 
theorem is proved and in an unpublished paper the same authors show that if n is an odd super 
perfect number, then ncr(n ) has at least 5 distinct prime factors, and that the number of distinct 
prime factors in n plus the number of distinct prime factors in a(n) is at least 7. 

94:13 (Carl Pomerance, via Doug E. Iannucci) For each positive integer n, is there a positive 
integer k such that a k (n)/n is an integer? 

Remark: a k (n ) = a(a k ~ 1 (n)), a°(n) = n. See the Cohen Sz te Riele paper mentioned above. 
94:15 (Paul Bateman) For a positive integer n let <f>„ denote the n-th cyclotomic polynomial, 

$ n (x) = I ^ 1 - x d Y {n/d) 

d\n 

Also let A n be the absolute value of the numerically largest coefficient in 4> n and let S n be the set 
of distinct integers occurring as coefficients in 4> n . If d and m are odd squarefree integers greater 
than 1 and if d, \ m, 

(a) is it always true that A d < A rn ? 
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(b) is it always true that S d C S rn ? 

Remark: Robert Harley, in email of 8 Feb 96, quotes Bateman to the effect that John 
Thompson originally asked a more specific question, namely, is it true that A(d) > 1 implies 
A(m) > 1 for all m which are multiples of d. If one could find a number m which (a) is a product 
of four distinct odd primes and (b) has the property that A(m) = 1, then there is a good chance 
that A(m/p ) > 1 for at least one of the four primes p dividing m. 

Harley continues: A few months ago I invented a new algorithm for computing cyclotomic 
polynomial coefficients quickly and applied it to search for such m. A search up to 1853733 found 
S m = { 0, ±1} when m is in 

{ 431985,757335,865365,1134915,1296885,1297815,1675365,1729335 }. 

For each p \ m, the coefficients of <1* (to jp) are again between -1 and 1. I initially expected to find 
counter examples, but now it looks like the answer is “yes” and I’m looking for a proof! There is 
a very interesting pattern in the factors of these m: 

431985 = 3 x 5 x 31 x 929 
757335 = 3 x 5 x 29 x 1741 
865365 = 3 x 5 x 31 x 1861 
1134915 = 3 x 5 x 29 x 2609 
1296885 = 3 x 5 x 31 x 2789 
1297815 = 3 x 5 x 31 x 2791 
1675365 = 3 x 5 x 61 x 1831 
1729335 = 3 x 5 x 31 x 3719 

In each case, a given factor plus 1 (or else minus 1) is divisible by all the smaller factors. This 
might lead to some insight for a proof. 

94:26 (Mike Zieve) For p prime, let 

f(x ) = x + bo + b ix + b 2 x 2 H-1- b n x n G ¥ p [b 0 , h,... , b n ][x] 

be a generic polynomial of degree n. Let g(x) G F p [b 0 ,... , b n \[x] be the p-th iterate of /, and 
write g(x) — x as a sum of monomials cx k bo°bi 1 ■ ■ ■ where c G F p \{0}, and k, a 0 , aq, ..., 
a n > 0. 

Conjecture: If k < p — 1, then 2ao + cq > p — 1. 

Remark: In April 1995 the proposer stated that Hendrik Lenstra had proved the conjecture 
and that the proof would appear in his, Mike Zieve’s, thesis. 
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Problems Proposed 16 & 19 December 1995 


95:01 (Joel Brenner and Richard Guy) Which primes occur as factors of numbers of the 
form 6 n + 7 n (or, more generally, of numbers satisfying a given second-order constant- 
coefficient linear recurrence)? Say something about the distribution of the residues of 
gn _|_ yn more generally, ...) modulo a prime. 

Remarks: 1. Concerning the first question, Basil Gordon and Carl Pomerance refer 
to a paper of Hasse; I think one they have in mind is 

Helmut Hasse, Uber die Dichte der Primzahlen p, fur die eine vorgegebene ganzrationale Zahl a#0 von 
gerader bzw. ungerader Ordnung mod. p ist, Math. Ann. 166 (1966) 19-23, MR 34 #5800, Revnum A16-52. 

Carl also writes out the following solution, based on Hasse’s argument. 

Theorem. The set S of primes p which divides some 7 n + 6 n has asymptotic density 
2/3 in the set of all primes. Among the density 1/3 of primes which never divide a 7 n + 6 n 
is a subset of density 1/4 (namely 3/4 of the primes which never divide a 7 n + 6 n ) which 
he in 1/4 of the residue classes mod 168 coprime to 168. 

Proof. Let £(p) be the order of 7/6 in the group (Z/pZ)* (assume p /42). Then 
7 n + 6 n = 0 (modp) (7/6) n = —1 (modp). Thus S = {p : £(p) is even}. Let 

T = {p : £(p) is odd}, so T is the set of primes which can never divide a 7 n + 6 n . Let 
Tk — {p £ T : 2 k || p — 1}. Then Ti = {p = 3 (mod 4) : (y) = 1}. Tj is a set of primes 
of density 1/4 defined via quadratic reciprocity by p lying in residue classes mod 168. By 
the Chebotarev density theorem, the density of T k is 4 -fe . Note that Ylk >l 4 _fe = 1/3- All 
it remains to note is that for each fixed A, the density of U k>KT k is at most 2 1 ~ K . since 
is contained in { p : p = 1 (mod 2 K ) }. Thus it is correct to add densities, since 
the density of the “tail,” namely U k>xTk, becomes small as K —* oo. 

2. Pieter Moree writes, THE reference for the first part of question 95:01 is: 

C. Ballot, Density of prime divisors of linear recurrences, Mem. of the AMS, 551, 1995. 

Pieter has two preprints which have a bearing on the question, a discussion of the 
history, and a good set of references. He continues: As to the second part of the question, 
here I would phrase the question as: on average how many residue classes (mod p) are 
covered by 6 n + 7 n (mod p|f So if 6 n + 7 n (mod p) assumes r(p) values, what is the 
average value of r(p)/p? I expect this to be some positive number. I expect that techniques 
as first explored by Hooley in his conditional proof (under GRH) of the Artin primitive 
root conjecture might yield the answer. 

95:02 (Peter Montgomery) Find a parametric solution for the diophantine equation xpiji + 
x 21/2 + xsys = 0. Are there quadratics ./i ■ f 2 , ./. 3 , 9i ; 92 ■, (js in indeterminates z 6 

such that “most” solutions of x\y\ + xpipi + xpy-s = 0 have the form Xj = f : j{z\...., zq), 
Dj = gj(zi..... zq). and such that fj9j identically zero? 

It will not do to let / 1 , / 2 , / 3 , hi, h 2 , h 3 be linear in Zi,.... z G and let ( 91 , 92 , 93 ) = 
(/ 1 , / 2 , /a) x (hi, h 2 , h 3 ), because the p 3 will typically be much larger than the fj. 

95:03 (J. C. Lagarias) Prove or disprove: every algebraic number held K has a finite 
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extension L whose ring of integers has a power basis. 

Remark: This is true for abelian number fields, where we can take L to be a cyclo- 
tomic field. 

95:04 (Neville Robbins) For i = 1,2 let F t (n) be the number of partitions of n into parts 
congruent to ±i (mod 5). Is it true that lim n ^ 00 = 1 — 

Solution: Basil Gordon answered this in the negative, proving that the limit is in 
fact (-s/5- l)/2. See 95:12. 

95:05 (Peter Borwein and David Boyd) Find a polynomial with coefficients in { 0, ±1 } 
and a 5th order zero of modulus strictly less than 1. More generally, show that such a 
polynomial can have a zero of arbitrarily high multiplicity strictly inside the unit circle. 
There is a connection to Lehmer’s problem on Mahler’s measure. 

95:06 (Paul Erdos, via John Selfridge) Let N be relatively prime to 6. Then N can be 
written as a sum of numbers of the form 2 a 3 b , no one dividing any other. 

1. Is there a function f(s) such that the least summand can be made greater than s 
for all N > f(s )? Thus, /( 1) = 1, /(2) = 5, /( 3) = 11; is it true that /(4) = 19, /(6) = 31, 
/(8) = 47, ...? 

2. Is it true that for every N the expression with maximum least summand is unique? 
E.g., for N — 115, the least summand must be at least 24, and the only expression with 
least summand 24 is 115 = 24 + 27 + 64. 

Solution to 1: Claudia Spiro-Silverman solves the first question; at a page and a 
half, the solution is a bit too long to be included here. It ends as follows: 

The function f(s) does, in fact, exist. And, the effectiveness of determining the 
value of the function /(s) as s increases is as good as the effectiveness of algorithms for 
determining 


max{ K : there is a power of 2 between 3 K — d and 3 K + d}, 

as a function of d. Thus, we were able to verify that /(3) = 11 by checking that the last 
solution of 3 k — 2 L = 1 is K = 2, L = 3. 

[Editor’s query: would Ellison’s result, 

\2 X — 3 y \ > 2 x e~ x/1 ° for x > 27 

(see UPINT F23) be of any use here?] 

95:07 (Melvin Knight) Let F m be the m-th Fibonacci number. 

1. If m = 5 n — 1, n = 1,2,..., and m! > m, then 5 n | (F m i — F m ), so lim^oo F m 
exists in the 5-adics; what is it? 

2. If m = 5 n , n = 1,2,..., and m' > m, then 5 n | (q^r- — ^), so lim n _ KX) ^ exists 
in the 5-adics; what is it? 
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Solution to 1: Peter Montgomery notes that F 2 n — F rn F m _ 1 — F 2 n _ l = 1 if m is 
odd, and 5 n | F rn if m = 5 n , whence F 2 l _ 1 = — 1 (mod 5 n ); it follows that the first limit 
is V^-T- 

Also with m = 5 n , Peter asserts that for every a there is a 5-adic integer K a such 

that 

- 3 + 0 («+ 

He notes that K -1 = Ai, A 3 = —3Ai, A 5 = 5Ai, and asks whether K a = «(—l)( a_1 )/ 2 Ai 
for all o. 

95:08 (Gary Walsh) Are there any positive integer solutions to x 2 (x 2 + 1) + y 2 {y 2 + 1) = 
z 2 {z 2 + 1 )? 

95:09 (Heng Huat Chan) It is well-known that Y2mn=- oo e ~’ K ^ m2+n ^ = v / 7 r/r 2 ( 3 / 4 ), 
where r(z) = / 0 °° e _t A _1 dz. It was recently shown that 

p -2Tr(m 2 +mn+n 2 ) _ _ \ _ 77 ^ 

^ (i2)V8vVfrir 2 (3/4)- 

For Q(m,n) = am 2 + bmn + cn 2 , define Vg(g) = Ylm n=-oo q < ^ ( ' m,n ^. For which Q and q 
can we explicitly evaluate Vg(g)? A less ambitious question would be, find an identity 
analogous to the given ones, involving ^ q m + mn + 2n . 

95:10 (Graeme Cohen and Herman te Riele) If 1 < a < b and <j)(a) = 0(6) = k > 3, 
can a be prime? 

Remarks: 1. If so, 6 has at least 12 distinct prime factors. 

2. 0(a) = 0(6) = has infinitely many solutions for k > 3, none for k = 1,2. 

3. a is prime if and only if (k — 1)0(6) = 6 + 1 — see UPINT B37 for some solutions 
of this equation. 

95:11 (Carl Pomerance) Suppose d\ < d 2 < ... and Yli=i di prime for all k. Show that 
there exists c > 0 such that d & > ck for all k. Must k~ x dk go to infinity? 

Remarks: 1. John Selfridge suggests a finite version (presumably, if d\ < d 2 < ■ ■ ■ < 
d n and Yli=i di is prime for 1 < k < n then bound n~ 1 d n from below). 

2. In email of 15 Feb 96, Carl writes: I was able to solve the first part (that d^ > ck) 
using the sieve. The second part is still open. Here is a related problem. Let gi < q 2 < ■ ■ ■ 
be any infinite sequence of primes. Is it possible for the second difference sequence for this 
sequence to have only finitely many values? Can the second difference sequence take values 
solely in the set 2,4? The second difference sequence is the sequence q n — 2g n _i + g n _ 2 . 

95:12 (Bruce Reznick) Suppose 0 < a\ < ... < a,j < m and 0 < 6i < ... < 6fc < n, 
and let p a (N) denote the number of partitions of N into parts congruent to a t (mod m) 
and similarly let Pb(N) denote the number of partitions of N into parts congruent to 6,; 
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(mod n). Under what circumstances does p a (N)/p b (N) converge to a finite non-zero limit 
as N —> oo, and how can the limit be expressed in terms of the parameters? 

Solution. Basil Gordon writes: It is convenient to replace the inequalities 0 < 
oi < ... < a,j < m and 0 < b\ < ... < bk < n by 0 < a\ < ... < a,j < m and 
0 < bi < ... < bk < n, respectively. An obvious necessary condition for p a (N)/p b (N ) to 
approach a non-zero limit as N —> oo is that gcd(ai, ... ,aj,m) = gcd(&i,..., bk, n) = 1. 
Assume henceforth that this holds. Then, as shown by Charles Brenner (J. Comb. Thy. A 
43 (1986) 303-319), we have 


p a (N)~a 0 N a 


where 02 = j/m, «i = j /2 — a i/ m , 


«o 


1 

2 y/n 


(1 r 2 « 2 \ 1/4 ~ ai/2 


TT £(flj/rra) 


and similarly for p b (N), with (3j in place of oij, j = 0,1,2 (here, T(s) is the gamma 
function). It follows that lirri/v^oc Pa(N) / p b (N) exists and is non-zero if and only if 0:2 = fh 
and cci = /?i, that is, j/m = k/n and 

j k 

j /2 — ^ di/rn = k/2 — ^ bi/n = 7 , say. 

i =1 i=l 

When these two conditions hold, the limit is 


^ = (n/m) 7 (27r)( fc 
Po 


ntiT {ai/m) 

nhm/n)' 


For example, taking m = n = 5, 01 = 2, 02 = 3, b\ = 1, 62 = 4, we have 02 = P 2 = 2/5 
and cci = /?i = 0. Hence the limit exists and equals 

u 0 = T(2/5)T(3/5) 

A) r(i/5)r(4/5)' 

We now apply the reflection formula T(.s)T(l — s) = 7 rcsc 7 rs to get 

ot 0 _ sin 7 r /5 1 \/5 — 1 

/3q sin27r/5 cos 7 t /5 2 


This solves 95:04. 

Gordon notes that simplifications are possible if one only wishes to solve the special 


case. 



95:13 (David Moulton) Find a nice expression for Yl p primed 2 - 

Remarks: 1. John Wolfskill suggests sending it to David Bailey (who has a pro¬ 
gram for finding linear relations among given constants), and Peter Montgomery suggests 
checking the Borweins’ dictionary of real numbers. 

2. In email of 12 Feb 96, Jeff Lagarias notes that 

iogc(s) = J2( p ~ s + \ p ~ 2a + \p~ 3s + • • •) 

p 

Now log£(2) gives the right first-order term, and log Q(2n) comes from 

inclusion-exclusion. 

3. Robert Harley independently gave the same formula for YhP~ 2 - Harley searched 
briefly for a linear relation with the familiar constants using lattice reduction, without 
success. He gives the sum to 500 decimals; the first 50 are 
0.45224742004106549850654336483224428476745164712962... 

95:14 (Richard Schroeppel) Consider a “computer” whose registers can store integers 
of arbitrary size. The computer can only do addition, subtraction, multiplication, and 
division with remainder; it can shift the read or write head to a different register, and has 
a finite state controller. 

Adi Shamir, Factoring numbers in 0(log n) arithmetic steps, Information Processing Letters 8 (1979) 28-31 

showed that one can factor N in polynomial time (« (log N) c operations) on such a 
computer. Is there an NP-complete problem which can be solved in a polynomial number 
of operations on such a computer? 

Remark: David Moulton suggested the partition problem; given non-negative integers 
fci,..., k n , can { 1,..., n } be partitioned into two sets S , T with Ylies ^* = This 

can be done if and only if f[[(±/ci ± ... ± h n ) — 0, taking the product over all choices of 
signs, and perhaps this product can be reduced by judicious use of identities to one which 
Schroeppel’s computer can handle. 

Solutions: 1. Ernest S. Croot III, in email of 3 January 96, writes; 

Here is a solution to [Rich Schroeppel’s] problem: 

The problem of choice to encode is the subset-sum problem. So, suppose that we want 
to decide whether the number N can be written as a subset sum chosen from the numbers 
oi, a 2 ,.. •, 0 ^. We consider the polynomial f(x) = (1 + x ai ) * (1 + x a2 ) * ... * (1 + x ak ). 
f(x) has the following two properties: (1) f(x) can be evaluated at x in polynomial time 
in the lengths of ai, a 2 ,..., a*,, and (2) the coefficient of x m in the expansion of fix) is the 
number of ways of representing m as a subset sum chosen from oi, a 2 , • • •, Ofe- 

Since the total number of subsets is 2 fe , the coefficients in the expansion of f(x) sum 
to 2 k . Clearly then, /( 2 k ) -r 2 Nk = c N + c N+1 2 k + c N+2 2 2k + ... + q2^- 7V ) fc , where a 
is the number of ways of representing i as a subset sum chosen from oi, a 2 ,..., a*,, and 
l = oi + o 2 + ... + afc (recall that # is integer division). Thus, N can be represented as 
a subset sum of oi, a 2 ,..., if and only if the remainder upon division of /( 2 k ) 4- 2 Nk 
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by 2 k is non-zero, and clearly this can be computed in polynomial time on the machine in 
question. 

2. Jeff Lagarias, in email of 18 Jan 96, adds this: 

Shamir was not aware that it had already been shown that register machines can solve 
any problem in PSPACE (which contains all of NP) in a polynomial number of steps, by 

J. Simon, On some central problems in computational complexity, Ph.D. thesis, Computer Science, Cor¬ 
nell 1975. 

It is also implicit in the paper: 

V. R. Pratt and L. J. Stockmeyer, A characterization of the power of vector machines, J. Comput. System 
Sci. 12 (1976) 198-221. 

The vector machine model is slightly different from the register machine model, but 
they are known to be of computationally equivalent power, according to David S. Johnson. 

The Shamir result shows the power of advertising. It was simple and easy to un¬ 
derstand, attractively presented, and got a great deal of publicity. The Simon result did 
not. 

95:15 (Doug Bowman) Let fo(x) = x, and, for x > n + 1, define f n+ i(x) by f n+ i(x) = 
fn+i f d (x) • there an efficient way to calculate f n {x)l for example, to find ./’iooo(lf)lO)? 

95:16 (Basil Gordon, via Gerry Myerson) For fixed n > 3, are there infinitely many 
nxn magic squares whose entries are distinct primes? What is wanted is an unconditional 
proof using, say, sieves, rather than one depending on unproved hypotheses such as, say, 
the prime /e-tuples conjecture. 

95:17 (Gerry Myerson) Let the matrix A be given by 


/ 2 

6 

14 

25 

45 \ 

21 

9 

18 

33 

11 ' 

24 

27 

28 

5 

8 

15 

43 

19 

3 

12 J 

V 30 

7 

13 

26 

16 j 


It is possible to write A — CjPj with ci,..., c r positive integers, P\. ... . P r permutation 
matrices, and r — 11. Is it possible with r = 10? The point is that there is such an 
expression with r = 10 for 2 A. 

95:18 (Martin LaBar, via Richard Guy) Is there a 3 x 3 magic square with distinct square 
entries? 

Remarks: This was problem 270, Coll. Math. J. 15 (1984) 69, reprinted as D15 in 
UPINT2. There is a discussion in the December 1995 Monthly, and John Robertson, 560 
Bair Road, Berwyn PA 19312, has a paper in Math. Mag. which may not have appeared 
yet. 

Robertson, in a letter to Richard Guy of 4 December 94, relates the problem to that of 
finding three Pythagorean triangles with the same area, the squares of whose hypotenuses 
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are in arithmetic progression, and to that of finding three points on an elliptic curve 
y 2 = x 3 — n 2 x which are the doubles of rational points and have their ^-coordinates in 
A.P., where n is a congruent number (see D27 in UPINT). 

Andrew Bremner finds six of the required eight magic sums in 

/15 2 20 2 60 2 \ 

36 2 48 2 25 2 
\52 2 39 2 0 2 ) 

Michael Schweitzer (eiffel@swissoft.li.provi.de) shows that any such square must have 
entries with at least 9 decimal digits. He gives the following specimen in which only one 
diagonal fails: 


127 2 

46 2 

58 

2 2 

113 2 

94 

74 2 

82 2 

97 


Bremner also gives a square of squares over Q (t), with only the principal diagonal 
failing. The entries are, 


f 2 ( 1249 - 2032f 2 + 2824f 4 - 4576f 6 - 392f 8 + 2624f 10 - 992f 12 + 128f 14 + 16f 16 ) 2 , 
4(—2 + 8 1 2 + f 4 ) 2 (-1 - 8 1 2 + 2t 4 ) 2 (5 - 2 1 2 + 2f 4 ) 2 (-7 + 4 1 2 + 2f 4 ) 2 , 

4f 2 (—1 - 8 1 2 + 2f 4 ) 2 (5 - 2t 2 + 2f 4 ) 4 (—7 + At 2 + 2f 4 ) 2 ; 

4t 2 (—1 - 8t 2 + 2t 4 ) 4 (5 - 2 1 2 + 2f 4 ) 2 (—7 + At 2 + 2f 4 ) 2 , 

t 2 (1201 - 2728 1 2 + 1168f 4 + 2384f 6 + 664t 8 - 2272t 10 + 832f 12 - 64t 14 + 16f 16 ) 2 , 

4(—1 - 8 1 2 + 2f 4 ) 2 (5 - 2 1 2 + 2t 4 ) 2 (—7 + At 2 + 2f 4 ) 2 (-2 - At 2 + 7f 4 ) 2 ; 

4(—1 - 8 1 2 + 2f 4 ) 2 (5 - 2 1 2 + 2f 4 ) 2 (—7 + At 2 + 2t 4 ) 2 (2 - 2 1 2 + 5f 4 ) 2 , 

4t 2 (—1 - 8t 2 + 2t 4 ) 2 (5 - 2 1 2 + 2f 4 ) 2 (—7 + 4t 2 + 2t 4 ) 4 , 

i 2 (—1151 + 3488f 2 + 1240t 4 - 5632t 6 + 3448t 8 - 256t 10 + 352f 12 - 256t 14 + 16t 16 ) 2 . 
It’s implicit in the work of Carmichael that there can be no 3 x 3 magic square with 
entries which are cubes or are fourth powers. 

95:19 (Graeme Cohen) Let be P the set of perfect numbers, M the set of multiperfect 
numbers (excluding 1), and H the set of harmonic numbers (excluding 1); 

P = { n = 2 }, M = { n > 1 ^ e Z }, H = { n > 1 : e Z }. Then M and H 

each strictly contain P, and no odd member of either set is known. Give a “sensible” set 
which contains M and/or H strictly and for which it appears to be difficult to find an odd 
member. The set 2Z will be accepted as an answer only if accompanied by a proof that it 
contains M and/or H. 
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